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ABSTRACT 

We present optical and near-IR identification and spectroscopy of the host galaxy of GRB 980613. 
The burst was apparently associated with the optically (restframe UV) brightest component of a system 
of at least five galaxies or galaxy fragments at a redshift of z = 1.0969. The component we identify as 
the host galaxy shows a moderately high unobscured star formation rate, SFR ~ 5 Mq yr _1 , but a high 
SFR per unit mass, indicative of a starburst. The image components show a broad range of (R — K) 
colors, with two of them being very red, possibly due to dust. Overall morphology of the system can be 
naturally interpreted as a strong tidal interaction of two or more galaxies, at a redshift where such events 
were much more common than now. Given the well established causal link between galaxy mergers and 
starbursts, we propose that this is a strong case for a GRB originating from a merger-induced starburst 
system. This supports the proposed link between GRBs and massive star formation. 

Subject headings: cosmology: miscellaneous — cosmology: observations — gamma rays: bursts 



1. INTRODUCTION 

Studies of the cosmic gamma-ray bursts (GRBs) are 
currently one of the most active and exciting fields of 
astrophysics. A deluge of important clues in the long- 
standing puzzle of GRBs began with the discovery of long- 
lived X-ray afterglows by the BeppoSAX satellite (Costa 
et at 1997). This was follo wed by the discovery of optical 
(van Paradijs et al. 1997) and radio afterglows (Frail et 
al. 1997), and then the first redshift measurement which 
demonstrated the cosmological nature of GRBs (Metzger 
et al. 1997). To date, about a dozen redshift measurements 
of GRBs have been obtained. Studies of afterglows con- 
firmed the syn chrotron shock model (e.g., Wijers, Rees & 
Meszaros 1997), and their physics now seems to be reason- 
ably well understood. For a recent review of observational 
results, see, e.g., Kulkarni et al. (2000). 

What is still not known is what are both the nature of 
the progenitors and the trigger that causes GRBs. The 
two currently popular models include mergers of neutron 
stars or other compact stellar remnants, leading to a cre- 
ation of a black hole, and a collapsar or hypernova model 
in which an explosion of a massive star produces a black 
hole remnant. In both cases, spin energy of a debris torus 
or of the black hole itself is used to power the GRB. While 
other models are still viable, there is now a growing evi- 
dence in favor of the collapsar/hypernova model, at least 
for the long-duration bursts detected by BeppoSAX. The 
most direct evidence favoring this was the probable detec- 
tion of sup ernova? associated with GRB 980326 (Bloom 



et al. 19991) and GRB 970228 ( [Reichart 1999| , |Galama ex 
al. 2000), as well as th e still controversial case of GRB 
980425 an d SN 1998bw ( |Galama et al. 1998| , [Kulkarni e\ 
al. 1998). All models involving massive stars and their 



remnants suggest that GRBs should be closely related to 
the massive star formation in galaxies. 



Study of GRB host galaxies can provide valuable clues 
which can constrain the models, and redshifts are nec- 
essary in order to derive the physical parameters of the 
GRBs itself, primarily the energy scale, as well as the 
hosts: their luminosities, star formation rates, morphol- 
ogy, locations of GRBs within them, etc. 

In this Letter we present deep imaging and spectro- 
scopic observations of the host galaxy of GRB 980613. 
GRB 980613 was detec t ed and localized by BeppoSAX 
QPiro, Costa et al. 1998| , [Piro et al. 1998| ). Following sev- 
eral unsuccessful attempts an optical transient (OT) was 
discovered by Hjorth et al. (1998). The detection of the 
host galaxy was reported by Djorgovski et al. (1998a). Its 
redshift determination was first reported by Djorgovski et 
al. (1999), and is described in more detail here, with ad- 
ditional data. 

2. OBSERVATIONS AND DATA REDUCTIONS 

In what follows we assum e the Galactic foreground ex- 
tinction of Ay — 0.27 mag ( Schlegel, Finkbeiner fc Davii 



1998[) , and use the Galactic extinction curve from Cardelli, 
Clayton, & Mathis (1998) with R = A v /E B -v =3.1. 

Our initial imaging date were obtained on the 
W. M. Keck Observatory 10-m telescope (Keck II) by 
Dr. H. Ebeling, on 16.30 June 1998 UT, in the R band us- 
ing the Low -Resolution Imaging Spectrometer (LRIS; Pkq 



et al. 1995| ); two 300 s exposures were obtained. We re- 



duced these images in the standard manner and confirmed 
that no cosmic rays affected the immediate area of the 
transient. After the transient had faded we reobserved the 
field in R band (29 November 1998 UT; 900s), I band (24 
March 1999 UT; 1000s), and K band (7 February 1999 
UT; 2040s). The final combined images in each filter is 
shown in fig. 1. The morphology of the system surround- 
ing the GRB is complex. We label 5 distinct components 



1 Partially based on the observations obtained at the W. M. Keck Observatory which is operated by the California Association for Research 
in Astronomy, a scientific partnership among California Institute of Technology, the University of California and the National Aeronautics and 
Space Administration. 
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2 The Redshift and Host Galaxy of GRB 980613 
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Fig. 1. — Late time images of the field of GRB 980613 with the Keck Telescopes in _R (left), I (middle) and K bands (right). Field size 
is 37" X 37", with N up and E to the left. The "x" marks the location of the OT. The offset star used for our spectrosopic observations is 
labeled "S" (the OT was 21.0"N and 3.9"E from the star), as is the Halpern et al. (1998) star "2" used for photometric zeropointing of our 
R and / band images. 



in fig. 2. Table 1 provides their magnitudes and offsets 
relative to the position of the OT. 

We referenced our R and / band magnitudes using star 2 
(Halpern et al. 1998; Diercks et al. 1998). For each compo- 
nent we used a circular aperature and performed an aper- 
ature correction using the curve-of-growth of star 2. The 
total magnitude of system A+B+C is R = 23.11 ± 0.05. 
We estimate the zeropoint systematic uncertainty due to 
an uncertain color correction at 0.15 mag. The 2-a up- 
per limit to a point source detection is R = 25.9 mag and 
/ = 24.9 mag. Our K band imaging was taken under pho- 
tometric conditions and we used SJ 9134 (Persson et al. 
1998) to obtain a zeropoint for the night. The K band 
aperature magnitudes of the 5 components are listed in 
Table 1. 

We registered the images from 16 June 1998 and 29 
November 1998 R band images using 7 unsaturated stellar 
objects common to both images within 30 arcsec from the 
host galaxy. We used the optimum-filter technique in the 
IRAF package CENTER to position the astrometric tie ob- 
jects. We determined the rotation, shift, and relative scale 
of the two images using GEOMAP and then registered the 
early-time image to the late-time image. The peak of the 
OT and the putative host were then estimated. Including 
the peak center errors and the registration uncertainty we 
find the OT was offset from the brightest R band peak (A) 
by 0.52"± 0.13" East and 0.83"± 0.14"S. This amounts to 
a radial angular projected offset of 0.98"± 0.14", and we 
interpret component A as the most likely host galaxy of 
the GRB. The / and K band images of the host were also 
registered to the late-time R band image. Fig. 2 shows a 
color composite 15.3" x 15.3"region around the GRB and 
its host using the late-time R, I, and K band images. 
The ellipse is a 3-<r error contour about the position of the 
GRB. 

Spectra of the host galaxy were obtained using LRIS 
on 17 December 1998, 14 January 1999, and 16 February 
1999 UT, all in good observing conditions. We used a 
1.5 arcsec wide long slit, always at a position angle close 
to the parallactic. On December 14, we used a 300 lines 
mm -1 grating, giving an effective instrumental resolution 
FWHM « 16Aand an approximate wavelength coverage 
3950-8950 A, and obtained 6 exposures totaling 8300 s. 



On January 14 and February 16, we used a 600 lines mm 
grating, giving an effective instrumental resolution FWHM 
w 8 Aand an approximate wavelength coverage 5970-8540 
and 5700-8270 A, respectively; 3 exposures of 1800 s were 
obtained on each night. The object was dithered on the 
spectrograph slit by several arcsec between the exposures. 
Exposures of an internal flat-field lamp and arc lamps 
were obtained at comparable telescope pointings imme- 
diately following the target observations. Exposures of 
standard stars from Oke & Gunn (1983) and Massey et 
al. (1998) were obtained and used to measure the instru- 
ment response curve. We estimate the net flux zero-point 
uncertainty, including the slit losses, to be about 10-20%. 




Fig. 2. — Color composite 15.3" X 15.3" region around the GRB 
and its host. We have used the late-time R, I, and K band images 
for the red, green, and blue channels. Before combining, we sub- 
tracted the mode-determined sky level and scaled each image by the 
r.m.s. noise of the background. The ellipse is a 3-cr error contour 
about the position of the GRB. The five distinct components A — E 
are labeled. The physical scale was is computed assuming all objects 
are in a redshift sheet at z = 1.096. As discussed in the text, the 
GRB appears within the light of component "A" , which we interpret 
as the host galaxy. 
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Fig. 3. — The weighted average spectrum of the host galaxy of GRB 980613 (component A in fig. 2), obtained at the Keck II telescope. 
Emission lines [O II] 3727 and [Ne III] 3869 are labeled; "ns" refers to noise spikes from a poor subtraction of strong night sky lines. 



Wavelength solutions were obtained from arc lamps in 
the standard manner, and then a second-order correction 
was determined from the wavelengths of isolated strong 
night sky lines, and applied to the wavelength solutions. 
This procedure largely eliminates systematic errors due to 
the instrument flexure, and is necessary in order to com- 
bine the data obtained during separate nights. The final 
wavelength calibrations have the r.m.s. ~ 0.2 — 0.5 A, as 
determined from the scatter of the night sky line centers. 
All spectra were then rebinned to a common wavelength 
scale with a sampling of 2.5 A using a Gaussian with a 
a = 2.5 A as the interpolating/weighting function. This 
is effectively a very conservative smoothing of the spec- 
trum, since it is smaller than the instrumental resolution. 
Individual spectra were extracted and combined using a 
statistical weighting based on the signal-to-noise ratio de- 
termined from the data themselves (rather than by the 



responds to 8.8 proper kpc in projection. 

The gamma-ray fluence from this burst w as (1.71 ± 
0.25) x 10~ 6 erg cm - 2 flWoods et at. 1998|). The cor- 



responding isotropic 7-ray energy is Ej^ so = 6.2 x 10 51 



erg. 



From the [O II] 3727 line flux, we derive the line lumi- 
nosity L3727 = 4.0 x 10 41 erg s _1 . Using the star forma- 
tion rate estimator from Kennicutt (1998), we derive the 
SFR w 5.6 Mq yr . From the UV continuum luminosity 
at \ res t = 2800A, following Madau, Pozzetti & Dickinson 
(1998), we derive SFR « 3.1 M© yr" 1 . The difference may 
be due to the internal reddening within the host galaxy, 
but we also note that the [O II] line estimator is more 
sensitive to the current or recent massive star formation 
than the UV continuum estimator. This is consistent with 
the presence of the [Ne III] 3 869 line, also seen i n spec 



expos ire time) . 

The hnal combined spectrum ot the galaxy is shown m 
fig. 3. A strong [O II] 3727 line emission is present, and 
a weaker [Ne III] 3869 line is also seen. The weighted 
mean redshift is z = 1.0969 ± 0.0002. The observed con- 
tinuum is moderately blue, with the spectrophotometric 



tra of some other GRB hosts (Bloom et at. 199S, Bloom. 
Djorgovski & Kulkarni 2001) which may be indicative of 
a recent, very massive star formation. We are completely 
insensitive to any fully obscured star formation compo- 
nent, so that these numbers represent lower limits. The 
only GRB host galaxy with a higher unob scured SFR mea- 



colorsl (B - V) w 0.45 and [V - R) w 0.35 mag. 

The corrected [U 11] 3727 line flux is (b.2h±0.1h) x 10 17 

erg cm s _1 Hz -1 , and its observed equivalent width is 
W\ — 125 ± 5 A, i.e., 60 ± 2.4 A in the restframe, which 
is on a high side for fiel d galaxies at compa rable redshifts, 
but not extraordinary ( [Hogg et at. 19981) . The [Ne III] 
3869 line has a flux of about 10% of the [O II] 3727 line. 
The continuum flux at \ bs = 5871 A, corresponding to 
Xrest — 2800 A, is F v = 0.84 fiJy, uncertain by a few 
percent (statistical). Additional spectrophotometric flux 
zero-point uncertainty is estimated to be ~ 10 — 20%. 

The continuum flux from our direct photometry in the / 
band, which corresponds roughly to the restframe B band, 
is F v — 1.15 /iJy, uncertain by ~ 10%. 

3. DISCUSSION 

We will assume a flat cosmology with Ho = 65 km s _1 
Mpc" 1 , fl M = 0.3, and A = 0.7. For z = 1.0969, the 
luminosity distance is 2.461 x 10 28 cm, and 1 arcsec cor- 



sured to date is the host of GRB 980703 (Djorgovski et at. 



1998) 



From the observed continuum flux in the restframe B 
band, we derive Mg ks —19.85 mag, which is about 1 mag 
fainter than the present-day galaxy. Considering that 
an average galaxy at this redshift may be ~ 0.5 — 1 mag 
brighter than today due to normal evolution effects, we 
conclude that this galaxy (A) is moderately underlumi- 
nous. The star formation rate per unit mass is thus fairly 
high, consistent with the large equivalent width of the [O 
II] line. We thus conclude that this galaxy is undergoing 
a starburst, albeit mild. 

However, the most interesting feature may be the mor- 
phology of the entire system ( A-E) , which is highly sugges- 
tive of an interaction or early stages of a merger of at least 
two galaxies, some of which may be partly obscured. Red- 
shifts of all 5 components are necessary to really test this 
interpretation; we note, however, that we detect a weak [O 
II] line emission coincident with the component C in our 
long-slit spectra, at the same redshift as the component A. 
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Another possibility is that we are seeing a chance superpo- 
sition of galaxies at very different redshifts, e.g., that the 
IR-bright components C and D may be unrelated to A (and 
perhaps also B and E), and if they are in the foreground, 
some gravitational lensing may be involved. However, we 
consi der this less likely than the interaction/merger hy- 
pothe sis. 



The rate of galaxy int eractions and merger s increases 
sharply with redshift (c/., Le Fevre et al. 2000 ), and find- 
ing a strongly interacting system at z ~ 1.1 is not surpris- 
ing. The projected separations of components A-E (see 
fig. 2) are typical of intergalactic separations where dark 
halos are overlapping and the tidal friction inevitably leads 
to a merger. 

The very red colors of the components C and D, (R — 
K) 4.4 and > 5.8, respectively, is naturally explained 
due to obscuration by dust, which is consistent with the 
merger hypothesis. Alternatively, they could be passively 
evolving ellipticals which formed at a very high redshift, 
which requires some fine-tuning of model parameters. 



It is now well established that mergers of gas-rich galax- 
ies lead to bursts of star formation. This may be the 
strongest case so far for a GRB from such a system. The 
only other published case to date, showing morphology 
consistent with a mild tidal interaction is t he host galaxy 
of GRB 990123 (e.g., |Bloom et al. 1999|, |Fruchter ' et al 



1999| . [Holland fc Hjorth 1999| ). This provides a further ev- 



idence in support of the connection of GRBs with massive 
star formation. 
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TABLE 1 

Photometry and Offsets of Components Near GRB 980613 



Comp. 


Offset (") 


R 


I 


K 




w.r.t. OT 


(mag) 


(mag) 


(mag) 


A 


0.52 W, 0.83 N 


23.81 ±0.06 


23.42 ±0.10 


21.65 ±0.22 


B 


1.58 E, 1.18 S 


24.95 ±0.17 


24.12 ±0.16 


> 22.3 


C 


0.25 E, 2.26 S 


24.45 ±0.10 


23.81 ±0.12 


20.08 ± 0.06 


D 


0.47 W, 3.84 S 


> 25.9 


> 24.9 


20.12 ±0.06 


E 


1.72 E, 3.44 S 


25.2 ±0.3 


24.66 ± 0.28 


> 22.3 



